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ABSTRACT 

O 

^ \ We present results from our analysis of the Infrared Space Observatory (ISO) 

^ \ data on the L 1641-N outflow region located in the Orion A molecular cloud. 

C^; The data were obtained with the array camera (ISOCAM) using two broad-band 

O \ filters, LW2 (6.7 /im), & LW3 (14 /im), and the narrow-band circular variable 

\ filter (CVF) which provides low resolution (R=A/5A ~ 40) spectra in the 5-17 /xm 

^ ' region. 

We detect a total of 34 sources in the 7.65' x 8.40' region covered by the 
^ \ broad-band filters. Four of these sources have no reported detection in previous 

\ studies of the region. The CVF data are available for only the central 3.2' x 3.2' 

portion of the total region, providing spectra for the brightest 7 of the sources in 
that region. 

We find that the source previously identified as the near-IR counter-part to the 
IRAS detected point-source (IRAS 05338-0624) is not the brightest source in the 
wavelength region of the IRAS 12 /im filter. We find instead that a nearby object 
(within the beam of IRAS and not detected at near-IR wavelengths) outshines 
all others sources in the area by a factor of ~2. We submit that this source is 
likely to be the IRAS detected point source. 

A comparison of the near-IR {J-H vs. H-K) and mid-IR {H-K vs. [6.7/im] — 
[14/im]) color-color plots shows that while at near-IR wavelengths only four (4) 
of the sources show evidence for emission above the values predicted from photo- 
spheric emission alone (hereafter referred to as excess emission), at least 85% of 
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all sources show evidence for excess emission at mid-IR wavelengths. This result 
supports similar conclusions from L-band surveys. 

The CVF spectra suggest a range of evolutionary status in the program stars 
ranging from embedded YSOs to young disks. When combined with optical and 
near-IR age estimates, these results show active current star-formation in the 
region that has been on-going for at least 2 Myr. 



1. Introduction 

Intensive observational and theoretical studies of young, optically obscured star-forming 
regions over the last four decades suggest that stars like our Sun begin their lives in the 
midst of cold and dense molecular clouds. During the earliest phases of their evolution, 
these young stars are embedded behind 50-200 magnitudes of visual extinction , and likely 
develop extended bipolar outflows and jets in the early phases of their evolution towards 
the main-sequence (for a review, see Lada & Lada 2002). The Orion giant molecular cloud 
complex is the closest example of such a star-forming region that offers both massive and 
low-mass stellar population. 

The L 1641-N region is located in the southern part of the Orion A molecular cloud. 
Attention was especially drawn towards L 1641-N by IRAS with the detection of a bright 
point source (IRAS 05338-0624). Subsequently, Pukui (1989) reported discovery of two 
well-separated outflow lobes surrounding the IRAS source. Further optical and near-infrared 
imaging stiidics revealed an enhancement of stellar density near the IRAS source at near- 
IR wavelengths (Chen et al. 1996; Hodapp et al. 1993), located "among the highest 
concentrations of HH [Herbig-Haro] objects known anywhere in the sky" (Reipurth et al. 
1998). 

The presence of the HH outflows implies that the L 1641-N region is one of the youngest 
star-forming environment in the Galaxy. Dynamic age estimates of the outflow sources 
presently observed in the Galaxy constrain outflow ages to less than 10^ years (Fukui 1989; 
Fukui et al. 1993). Studying these young, optically obscured environments at mid-infrared 
wavelengths offers several advantages: (i) stars in the earliest phases of evolution are heavily 
obscured and are frequently undetected at even near-IR wavelengths, (ii) With imaging 
data alone, it is difficult to distinguish between actual photospheric emission and reprocessed 
emission (reflected hght). Thus, near-IR reflection nebulae are often confused with actual 
sources, (iii) At mid-IR wavelengths, both the stellar photosphere and the colder star-forming 
environment are visible. 
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However, the mid-IR region is difficult to observe from the ground because of the pres- 
ence of high background from the warm Earth atmosphere. The Infrared Space Observatory 
(ISO, Kessler et al. 1996) provided the first space-based observatory with the abihty to 
carry out sensitive observations at mid-IR wavelengths. The ISO observed spectra of the 
embedded Young Stellar Objects (YSOs) revealed, for the first time, the presence of silicate, 
CO, water (vapor and ice) absorption features in the photospheric continuum of the stars 
(see for example, Gibb et al. 2000). These features are thought to arise in the colder 
star-forming environment in which the stars are forming. Thus, the mid-IR region provides 
information both on the star itself and about the physical conditions in which they form. 

In this contribution we provide an analysis of ISO CAMera (ISOCAM, Cesarsky et al. 
1996) mid-infrared broad-band and low-resolution Circular Variable Filter (CVF) observa- 
tions of the L 1641-N region. Section 2 describes these observations in detail. Section 3 
describes our source detection, photometry and cross-correlation with other optical and 
near-IR sources. Wc present and discuss our results in Section 4. Finally, the conclusions 
are summarized in Section 5. 



2. Observations and data reduction 

All data for this contribution were obtained from the ISO Archive (Salama et al. 2002). 
The observing parameters and principle investigators for the observations used in this con- 
tribution are summarized in Table 1. Column 1 in Tabic 1 lists the unique ISO identifier, 
the TDT number, column 2 lists the ISO Astronomical Observing Template (AOT) used to 
obtain the data, columns 3 and 4 identify the filters and the associated wavelengths. For the 
LW2 and the LW3 broad-band filters, the central wavelength is listed, and for the Circular 
Variable Filter (CVF) the wavelength coverage is listed. Column 5 lists the total exposure 
time and column 6 lists the principle investigator (PI) for the observing proposal. 

The data were reduced using the CAM Interactive Analysis (CIA, Ott et al. 2001) 
software. The raw data files (CISP & IIPH files) were obtained from the ISO Archive and 
reprocessed locally. The data reduction comprised both "standard" CIA and customized 
reduction modules. The reduction procedure - in order of steps performed - is listed below. 

Dark Current Correction Over the lifetime of ISO, the ISOCAM dark frames showed cycli- 
cal variations in the observed levels of dark current. The periodicity, level and causes of 
this ISOCAM dark behaviour are described in Blommaert et al. (2001). The CIA routine 
'corr_dark' was used with the 'VILSPA' option to subtract the dark current from each in- 
dividual frame. The 'VILSPA' dark model accounts for all known ISOCAM dark current 
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variability modes. 

Cosmic Ray Removal The CIA routine 'deglitch' with the multi-resolution median filtering 
option was used to provide an initial assessment of the cosmic ray hits. This automatic 
procedure identified most cosmic ray events; however, a small (< 2%) of the pixels were 
incorrectly identified as cosmic ray hits. Thus, the results from the automated procedure 
were augmented by visual inspection of the data cube. 

Non-linearity corrections The ISOCAM detectors show the so-called "transient" behaviour 
- a non-linear response to input flux. The theory of detector transients is described by, for 
example, Haegel et al. (2001). The ISOCAM transients are reviewed by Coulais et al. 
(2002). Wc used the 'Fouks-Shubert' removal within CIA to correct for transients. This 
correction did not remove all transient artifacts in the broad-band rasters (where the effects 
are the largest). We simply masked out the uncorrected data and replaced these values by 
redundant values from the overlapping raster observation. 

Flat-field correction A single flat fleld for the entire array was constructed by median- 
combining all raster segments. This option is preferable to using a library flat-fleld, since a 
custom corrected flat-field is not aflected by the "jitter" in the ISOCAM optics, which moves 
the optical flat-field pattern on the ISOCAM array. 

The CIA routines 'raster _scan' was used to assemble the final mosaic of the broad-band 
data. 



3. Source detection and photometry 

All sources were identified using the DAOPHOT (Stetson 1987) find algorithm as 
implemented within the Interactive Data Language (IDL) environment. Visual inspection of 
the detected sources revealed that most DAOPHOT detections were either bright knots in 
the nebulosity or ISOCAM artifacts. The final sample of stars is, hence, drawn by confirming 
each detection with visual inspection or by comparison with the 2MASS catalog. Given the 
unreliable nature of source finding algorithm used - complicated by varying nebulosity and 
undersampled images, and uncertain PSF - there may be additional sources but these cannot 
be confidently seperated from non-sources. Hence, only those which are verified by one of 
the above methods are included in our final list. 

We used a 3 pixel by 3 pixel square aperture for photometry on all sources. We use the 
ISOCAM PSF library (available within CIA) to correct for flux outside of our adopted square 
aperture. This correction is 12% and 15% for LW2 and LW3 filters, respectively. We use a 
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global background sky estimate for all sources. To obtain the sky /background per pixel, we 
determine the median value in all regions that show no extended emission. Similar apertures 
were also used to extract the CVF spectra of the sources, where CVF data is available. No 
background correction is applied to the CVF spectra. 

We used the CIA's internal flux cahbration routine 'adutomjansky' to convert source 

counts to physical flux units. This calibration is based on standard stars and is discussed by 
Blommaert et al. (2001) and gives the following values for the LW2 and the LW3 filters: 

LW2: 1 mJy = 2.3303 counts 
LW3: 1 mJy = 1.9733 counts 

We refined the initial astrometry estimates for the sources (based on ISO's pointing 
information) by comparing the images to 2MASS images and using the 2MASS coordinates, 
and by applying the optical distortion correction given by Ali et al. (2003). Given the low 
number and uncrowded (source PSFs are well-separated) field, we cross-identified our ISO 
sources visually with 2MASS catalog and with the near-IR study of Chen et al. (1993). 

The final photometry and astrometry measurements for all sources discussed in this con- 
tribution is hsted in Table 2. Column 1 hsts our assigned identification (ID) number. Col- 
umn 2 list Chen et al. (1993) designations for the sources. Wc list multiple cross-identified 
sources when a single unambiguous determination could not be made. Columns 3 & 4 list 
the 2MASS coordinates for the source. Sources with ID numbers 10, 18, 19, & 20 have no 
identified 2MASS counterpart. For these sources the ISO coordinates are listed. To ensure 
uniformity, we computed coordinate offsets between ISO and 2MASS positions. For sources 
without 2MASS detections the hsted ISO coordinates are corrected for the aforementioned 
positional offsets. The magnitude of the positional offset is 0.24 and 7.7 arc-seconds in right- 
ascension and declination, respectively. Columns 5 and 6 list the ISO flux in milli-Jansky. 
The errors are computed from the standard deviation observed among the multiple mea- 
surements of the object. For eight of the sources, only one measurement is available. We 
assigned 20% errors to these sources, which is the average error percentage for all sources 
with multiple measurements. These stars are identified in Table 2. Four objects were dis- 
cernible only in the final mosaicked image of the region. For these objects, we used the 
XPHOT photometry package available within CIA (Ott et al. 2001) to obtain aperture 
photometry. These sources are also identified in Table 2. Columns 7-9 list the near-infrared 
flux measurements from 2MASS. 

We compared our photometry with that made available by the ISOCAM pipe-line (OLF 
version 10.0, Blommaert et al. 2003). We consider the OLP photometry to be less rehable 
in that several of the artifacts {e.g. residual images from transients. Section 2) are visible 
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in the OLP reduced data. These were removed from our reduction. Our comparison shows 
that mean differences between our photometry and OLP photometry to be 18% and 34% for 
the LW2 and the LW3 filters, respectively. 

4. Results and discussion 

We detect a total of 34 sources in the region covered by the broad-band filters. The 
spatial coverage and sensitivity of the narrow-band filter data allowed us to obtain low- 
resolution spectra of only 7 of these sources. All but 4 of the sources are previously reported 
by the Two Micron All Sky Survey (2MASS) and/or the near-IR survey of Chen et al. 
(1993). Two of the sources (^^10 and #18, in our numbering scheme) are only detected at 
mid-IR wavelengths: by ISOCAM and by Chen et al. (1993) in their M-band images. These 
results are further discussed below. 



4.1. The counterpart to IRAS 05338-0624 

Chen et al. (1993) suggest source #18, in our numbering scheme, as the likely near- 
IR counterpart to the IRAS detected point source based on their M-band imaging survey. 
Source #18 is within 2 arc-seconds of the IRAS published position and is the brightest source 
at M-band in the Chen et al. (1993) survey. Stanke, McCaughrean, & Zinnecker (1998) 
also detected a bright (0.4 Jy) point source at lO/xm coincident with the location of source 
#18. However, Stanke, McCaughrean, & Zinnecker (1998) argue that IRAS 05338-0624 
cannot be unambiguously associated with source #18 because a 2.7-mm dust continuum 
source (Chen et al. 1993, 1995) is also coincident with the position uncertainty ellipse of 
the IRAS source (Figure 2, also see Figure 4 of Stanke et al. . 1997). 

In the ISO image, source #10 outshines source #18 by approximately a factor of two. 
Source #10 is outside the field-of-view of the Chen et al. (1993) M-band survey and the 
Stanke, McCaughrean, & Zinnecker (1998) 10//m image, and is ~1 arc-minute away from 
the IRAS source position. There are no other comparably bright sources detected in the 
ISO images. The narrow-band CVF spectra (Figure 3) include the wavelength range of the 
shortest IRAS filter and are useful for obtaining a more appropriate fiux comparison with 
IRAS measurements. We used the ISO Spectral Analysis Tool (ISAP) and the IRAS filter 
transmission curves to simulate IRAS photometry for sources #18 and #10 using the CVF 
spectra. We find source fluxes of 284 and 481 mJy (30% errors) for sources #18, and #10, 
respectively. The IRAS published photometry for IRAS 05338-0624 is 481 mJy. 
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Despite the excellent agreement between the fluxes of IRAS 05338-0624 and source #10, 
we conclude that source #18 remains the most likely counter part to IRAS 05338-0624 be- 
cause it is the brightest source that is consistent with the position uncertainty ellipse of the 
IRAS source. Wc rule out source #10 because for it to be the IRAS counterpoint implies a 
positional error of ~1 arc-minute in the IRAS position. We consider this to be unlikely given 
that IRAS positions have been shown to be reliable to better than 7 arc-seconds (1-sigma) 
based on statistical comparison between IRAS and the Smithsonian Astrophysical Obser- 
vatory (SAO) catalog point sources (cf. Section VII of the IRAS Explanatory Supplement, 
Beichman et al. 1988). The 2.7-mm dust continuum source lacks a bright mid-IR counter 
part in the ISO image implying either that it is not a strong emitter at those wavelengths 
or that it has dimmed significantly in the decade between IRAS and ISO measurements. 

However, we also consider it unhkely that IRAS would not have spatially resolved and 
detected source #10, which is the brightest source in the ISO image and about 1 IRAS 
beam- width away from IRAS 05338-0624 (source #18). If wc rule out positional error in 
IRAS coordinates, then this implies that source #10 must have brightened by about at least 
200 mJy (using point source detection limit of 0.3 Jy for IRAS) between the ISO and IRAS 
epochs. 

We also note that the proposed IRAS counterpart, source #18, shows a broader PSF 
than other point sources in the area. This is likely due to the presence of extended emission 
or emission knots which are unresolved by ISOCAM and seen in reflected light in the K-band 
as noted by Chen et al. (1993). 

4.2. The nature of point-like sources 

Figure 4 shows the J-H vs. H-K, color-color relationship for the 30 sources for which 
near-IR photometry is available. In constructing this relationship, we used the photometry 
from 2MASS database. The solid lines in the figure show the colors of ordinary dwarf and 
giant stars from the compilation of Bessel & Brett (1988). The dotted lines in Figure 4 show 
the effect of extinction on the colors of the reddest giant and dwarf stars. The distribution 
of the L1641-N stars on Figure 4 is consistent with those from deeply embedded young 
stellar clusters. And, four stars, #6, #11, #17, and #21 show evidence for the so-called 
"excess' emission - colors that are redder than predicted by simple foreground extinction. 
This "excess" emission is though to arise from warm circumstellar material (see e.g. Lada 
& Lada 2003). The emission itself is from warm dust in the circumstellar material. 

We investigated a similar color-color relationship by using the ISO photometry for which 
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all 4 photometry measurements are available (24 stars). The results are shown in Figure 5: 
J-K vs. [6.7 fim] — [14/im]. As in Figure 4, the dotted line shown in Figure 5 shows the 
effect of foreground extinction, but in this case, for a star with intrinsic [6.7/im] — [14/im] 
color of mag. At mid-IR wavelengths only a handful of estimates are available for the 
intrinsic colors of ordinary or standard stars. These estimates show intrinsic values that 
are either or close to mag. Thus, we adopted mag as the typical [6.7/xm] — [14/im] 
color. While 4 stars show evidence for excess emission in Figure 4, all but four, sources #5, 
#31, #22, and #34 show evidence for excess emission at mid-IR wavelengths. This can 
be explained simply as warm dust emits more at longer than near-IR wavelengths; thus it 
is easier to detect excess emission at mid-IR wavelengths than near-IR color-color plot. A 
similar argument was presented for the Chamaclcon I population by Nordh et al. (1996) 
and Persi et al. (2000). Adding source #10 and #18 in this sample, we conclude that at 
least 85% of all sources detected in the ISO LW2 filter show evidence of excess emission at 
mid-IR wavelengths. 

Given the presence of excess emission in most of the observed sources, we conclude 
that most sources are most likely physically associated. A conclusion also reached by Chen 
et al. (1993) in their analysis, who labeled it as a "stellar density enhancement". The ISO 
narrow-band CVF spectra for the seven sources detected in the CVF filter shed further light 
on the nature and evolutionary status of the stars in this grouping. These ISO spectra are 
shown in Figure 3. Each panel of Figure 3 shows the spectrum from the source identified in 
the left-hand-side of the panel. Some of the prominent features associated with young stars 
are also labeled. 

Five of the 7 detected sources, #11, #13, #10, #18, and #19, show the broad feature 
at 9.5 //m, which is characteristic of sihcate absorption. Sources #10, and #19 further show 
features seen in the very youngest proto-stars such as W 33A (Gibb et al. 2000). These 
features are labeled in Figure 3 and have been attributed to ices (CO2, H2O, etc.) in the 
environment of the young embedded stars (Gibb et al. 2000). Source #14 shows marginal 
detections at wavelengths consistent with H2 emission. Sources #11 and #15, which show 
evidence for excess emission at near-IR wavelengths (see Figure 4), also show evidence for 
an emission feature at wavelengths consistent with the PAH emission at 11.2 yum. We rule 
out silicate emission because its peak is near 10 jim.. PAH emission has been observed in 
some T-Tauri stars and roughly 50% of the Her big Ae/Be stars in the sample of Meeus et al. 
(2001). PAH emission can signify the presence of tiny grains and (likely) UV radiation to 
excite the PAH. This is consistent with the observed location of sources #11 and #15 on 
Figure 4. Source #19, the faintest object detected in the CVF spectra, has relatively poor 
signal-to-noise ratio and most narrow features seen in its spectra are likely instrumental 
artifacts. The only significant detection in the spectra of source #19 is the broad silicate 
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feature. 

Evans et al. (2003) provide a crude evolutionary sequence for the youngest stars 
based on spectral features in the mid-infrared wavelengths (see their Figure 8) . We use this 
sequence to interpret our CVF spectra. Sources #10, #18, and possibly #19 are, thus, in 
the earliest phase of evolution, labeled as 'embedded YSOs' by Evans et al. (2003). The ice 
features in the spectra of sources #10, #18 are seen in the earliest phases when the proto-star 
is deeply embedded in a circumstellar envelope. As the circumstellar environment matures, 
the ices give way to molecular and atomic features during the 'Emerging YSO' phase. Source 
#14 in our survey is consistent with such an object. Sources #11, #15 show PAH emission 
but otherwise a featureless spectrum. These are consistent with the 'Young disk' phase of 
Evans et al. (2003). The evolutionary status of source #13 is unclear. The presence of 
sihcate absorption is indicative of the embedded YSO phase but the lack of features from 
other ices is puzzling. Perhaps this is an intermediate object between the Embedded and 
Emerging YSO phase. 

In our interpretation of the spectra the ages of stars are of the order ~fewxlO^ years for 
the embedded YSOs to ~fewx 10^ years for stars with 'Young disk'. These ages are consistent 
with Pukui (1989) based on their dynamical estimates using the outflows. However, Chen 
et al. (1993) estimate the cluster population to be 1-2 Myr based on 3 stars detected by 
optical and near-lR surveys. This inconsistency is likely just the difference between the 
optically detected and the mid-lR detected population of L 1641-N. The combination of 
the two estimates establish this region as one of the youngest regions in Orion with active 
star-formation and age spread of order 2 Myr among the stars. 

5. Summciry 

Our analysis of the ISOCAM broad-band and CVF data produced the following main 
results: 

We detect a total of 34 sources in the region covered by the broad-band filters. The CVF 
spectra are available for 7 of these sources. Two sources (#10 and #18) are only detected 
at mid-IR wavelengths. 

We find that the source previously identified as the counter-part to the IRAS detected 
point-source in the region is not the brightest source at the IRAS 12 urn filter wavelengths. 
We propose instead that the IRAS source is likely to be combination of sources with #10 
(in our listing) being the dominant one. This source outshines all others by roughly a factor 
of 2. 
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However, source #18, previously associated with the IRAS source, shows a larger PSF 
than other point sources in the area. This is likely due to the presence of extended emission 
or emission knots which are unresolved by ISOCAM. We conclude that source #18 remains 
the most plausible source of the outflows as suggested by Chen et al. (1993). 

A comparison of the near-IR {J-H vs. H-K) and mid-IR [H-K vs. [6.7/im] — [14//m]) 
color-color plots shows that while only two of the sources show evidence for "excess" emission, 
at least 85% of all sources show evidence for excess emission at mid-lR wavelengths. For 
two of these sources. CVF spectra arc available and show emission features from PAH. The 
presence of excess emission in most sources supports the assertion by Chen et al. (1993) 
that these objects arc physically associated, or at least co-located. 

The narrow-band CVF spectra of the 7 sources detected in the CVF show a range of 
evolutionary status for the stars, ranging from the earliest 'Embedded YSOs' to the 'Young 
disks'. When combined with optical and near-IR age estimates, these results show active 
current star-formation in the region that has been on-going for 2 Myr. 

The observed absorption features in the spectra of the sources deserve further study in 
the form of detailed modeling. Future work will include full radiative transfer modeling to 
interpret the physical conditions responsible for the observed spectral features. We plan to 
provide these results in a future contribution. 

We thank the anonymous referee for helping us improve this manuscript. A.N-C. re- 
search is partially supported by NASA through a contract with Jet Propulsion Laboratory, 
California Institute of Technology, and ADP Grant NRAOOOl-ADP-096. 
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Table 2. Photometry of L 1641-N sources. 
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Table 2 — Continued 
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^Assigned 20% error. See text for more detail. 
'^Photometry obtained from the final mosaic 
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Fig. 1. — The final reduced ISOCAM mosaic of the L 1641-N region. The image has been 
rotated such that North is to the top and East is to the left. The image size is 7.65 by 8.40 
arc-minutes. 



-17- 




Fig. 2. — Same as Figure 1 except with the identification labels for the sources identified. 
The region marked with the solid line shows the coverage of the CVF data. 
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Fig. 3. — The low-resolution spectra of 7 sources scanned with the CVF. Each object is 
identified in the left-hand side of the panel. The prominent spectral features are identified. 
These spectra are consistent with those seen for the youngest proto-stars in star-forming 
regions. 
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Fig. 4. — The J-Hvs. H-K color-color relationship for all ISOCAM sources also detected 
by 2MASS. The solid lines show the location of normal dwarfs and giants on this diagram. 
The dotted lines show the reddening vectors. Most of the L 1641-N sources appear reddened 
similar to other star-forming regions. Four sources show evidence for "excess", presumably 
circumstellar, emission. 
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Fig. 5. — The J-Kvs. [6.7//m] — [14/im] color-color relationship of all sources detected in all 
4 filters. The lines have the same meaning here as in Figure 4. All but four source seen 
show evidence for excess emission. Thus, at mid-IR wavelengths, at least 85% of the sources 
show excess emission, consistent with Persi et al. (2000) who find 75% of the sources in the 
Chamaeleon I cloud show excess emission at mid-IR wavelengths. 



